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ABSTRACT 

By means of high-resolution N-body simulations we revisited the dissipationless 
collapse scenario for galaxy formation. We considered both single-component collapses 
and collapses of a cold stellar distribution in a live dark matter halo. Single-component 
collapses lead to stellar systems whose projected profiles are fitted very well by the Sersic 
i?^/™ law with 3.6 < m < 8. The stellar end-products of collapses in a dark matter 
halo are still well described by the i?^/™ law, but with 1.9 < m < 12, where the lowest 
m values are obtained when the halo is dominant. In all the explored cases the profiles at 
small radii deviate from their global best-fit i?^/™ model, being significantly flatter. The 
break-radius values are comparable with those measured in 'core' elliptical galaxies, and 
are directly related to the coldness of the initial conditions. The dissipationless collapse 
of initially cold stellar distributions in pre-existing dark matter haloes may thus have 
a role in determining the observed weak homology of elliptical galaxies. 
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1 INTRODUCTION 

^ It is a well established fact that the end-products of dissipa- 
, tionless collapses reproduce several structural and dynamical 

■ properties of elliptical galaxies. For example, the pioneer- 
, ing work of van Albada (1982, vA82) showed that the end- 

■ products of cold collapses have projected density profiles well 
' described by the 7?^^* de Vaucouleurs (1948) law, radially de- 
creasing line-of-sight velocity dispersion profiles, and radially 
increasing velocity anisotropy, in agreement with observa- 
tions of elliptical galaxies (see also May & van Albada 1984; 
McGlynn 1984). More recently the dissipationless collapse 
has been studied in greater detail thanks to the advances in 
N-body simulations (e.g., see Aguilar & Merritt 1990; Lon- 
drillo, Messina & Stiavelh 1991; Udry 1993; Trenti, Bertin 
& van Albada 2004). These studies show that a smooth fi- 
nal density distribution with R^^^ projected mass profile is 
produced when the initial distribution is cold, extended, and 
clumpy in phase-space. From the astrophysical point of view 
the dissipationless collapse (Eggen, Lynden-Bell & Sandage 
1962) was introduced to describe a complex physical scenario, 
in which the gas cooling time of the forming galaxy is shorter 
than its dynamical (free-fall) time, stars form 'in flight', and 
the subsequent dynamical evolution is a dissipationless col- 
lapse. 

From the observational point of view, elliptical galaxies 
are characterized by a systematic luminosity dependence of 



their surface brightness profile, the so-called weak homology. 
For example, high quality photometry of elliptical galaxies 
indicates that the R^^"^ Sersic (1968) law (see equation 9 
below), with best-fit index m ranging from m ~ 2 for faint 
ellipticals up to m ~ 10 for the brightest ones, often repre- 
sents their surface brightness profiles better than the de Vau- 
couleurs law (the m = 4 case; Caon, Capaccioli & D'Onofrio 
1993; Prugniel & Simien 1997; Bertin, Ciotti & Del Principe 
2002; Graham & Guzman 2003). In addition, Hubble Space 
Telescope observations probing the luminosity profiles of sev- 
eral ellipticals down to sub-arcsec resolution (Ferrarese et 
al. 1994; Lauer et al. 1995) revealed that at these small radii 
some profiles fiatten (core galaxies), others are characterized 
by steep cusps (power-law galaxies). While the surface bright- 
ness profiles of power-law galaxies are well fitted by the R^^"^ 
law down to centre, within the break radius Rh the surface 
brightness profile of core galaxies stay well below their global 
best-fit Sersic model (Graham et al. 2003; Trujillo et al. 2004; 
Ferrarese et al. 2006). In general, core profiles are common 
among bright ellipticals, while faint systems tend to have 
power-law cusps; remarkably, other galaxy global properties 
are related to the presence of the core (e.g. Pellegrini 1999, 
2005). 

The explanation of the observed weak homology of el- 
liptical galaxies is important to understand galaxy forma- 
tion. For example, the presence of a core is usually inter- 
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preted as the signature of the merging of supermassive black 
holes, consequence of the merging of the parent galaxies (e.g. 
Makino & Ebisuzaki 1996; Faber et al. 1997; Milosavljevic 
& Merritt 2001), while in N-body simulations of repeated 
equal-mass dissipationless ('dry') mergers the best-fit m of 
the end-products increases with their mass (Nipoti, Londrillo 
& Ciotti 2003a, hereafter NLC03). However, Ciotti & van Al- 
bada (2001) and NLC03 showed that repeated dissipationless 
merging fails to reproduce the Faber & Jackson (1976), the 
Kormendy (1977), and the Mbh-o" (Gebhardt et al. 2000, 
Ferrarese & Merritt 2000) relations, and also that a substan- 
tial number of head-on minor mergings make m decrease, 
bringing the end-products out of the edge-on Fundamental 
Plane. These results, together with other astrophysical ev- 
idences based on stellar population properties such as the 
Mgj-cr relation (e.g. Burstein et al. 1988; Bender, Burstein 
& Faber 1993), indicate that dry mergings cannot have had a 
major role in the formation of elliptical galaxies, and gaseous 
dissipation is needed (see, e.g., Robertson et al. 2005; Naab 
et al. 2005). 

In alternative (or as a complement) to the merging sce- 
nario, it is then of great theoretical interest to explore if (and 
if yes under what conditions and to what extent) the dissipa- 
tionless collapse of the stellar population produced by a fast 
episode of gaseous dissipation and the consequent burst of 
star formation is able to reproduce end-products with pro- 
jected density profiles well described by the Sersic law. In 
particular, given that according to the current cosmological 
picture galaxies form at peaks of the cold dark matter dis- 
tribution (e.g. White & Rees 1978; White & Frenk 1991), 
it is natural to investigate dissipationless collapse in two- 
component systems. In this work we study this process using 
high-resolution two-component N-body simulations, in which 
the collapse of the stellar distribution and the response of the 
dark matter halo are followed in detail. As will be shown, 
dissipationless collapses in pre-existing dark matter haloes 
are indeed able to reproduce surprisingly well the observed 
weak homology of elliptical galaxies and the fiat inner sur- 
face brightness profiles of core ellipticals arise naturally from 
dissipationless collapse, with i?b determined by the coldness 
of the initial conditions. 

This paper is organized as follows. The numerical simu- 
lations are described in Section 2. The results are presented 
in Sections 3. Section 4 summarizes. 



2 NUMERICAL SIMULATIONS 
2.1 Initial conditions 

We consider two classes of simulations. In the first we follow 
the virialization of a cold, single-component density distri- 
bution. In the second the initial conditions represent a cold 
component (stars) deemed to collapse in a nearly-virialized 
dark matter halo. 

The initial conditions consist of the stellar (p,) and the 
halo (ph) density distributions, and the corresponding virial 
ratios /3, and Ph measure the 'coldness' of the distributions. 
For example, the stellar virial ratio is defined by 



2K, 



ii7, +W^,J 



(1) 



(2) 



are the kinetic energy and the self-gravity of the stellar dis- 
tribution, 



W.h = - / < X , V$h > p.d X 



(3) 



is the interaction energy between stars and dark matter, "I>* 
and $h are the gravitational potentials of stars and dark 
matter, respectively, and <, > is the standard inner product. 
Note that in one-component simulations /3t = 2Kt/\U,\. All 
simulations are evolved up to virialization, which is reached 
typically within 40fdyn, where 



tdy 



GM, 



5/2 



(2|£tot|)3/2 



(4) 



is the dynamical time of the system, and Aftot and -Etot are 
its total mass and energy, respectively. 



2.1.1 Stellar component 

The initial configuration of the stellar component is obtained 
by introducing inhomogeneity on a smooth, spherically sym- 
metric density distribution p*. In particular, we adopt 



p.(r) 



3M*r: 



47r(r2 + r|)5/2 
(Plummer 1911), and 

3-7 M,r* 



p.(r) = 



47r rT(r, + r)*-T 



(5) 



(6) 



with 7 = or 7 = 1 (Dehnen 1993; Tremaine et al. 1994), 
where M, is the total mass and r* the characteristic radius. 
The particles are spatially distributed according to equations 
@ or JSJl and then randomly shifted in position (up to rt/5 
in modulus), so that the distribution results inhomogeneous. 

The velocities of the particles are first randomly ex- 
tracted from an isotropic Gaussian distribution with a given 
value of the variance, and we measure the corresponding 
virial ratio f3o- The required value of the stellar virial ratio 
0.002 ^ /3* ^ 0.2 is then obtained by scaling the velocity 
of each particle by ^/3« //3o . In some case the stellar distri- 
bution is characterized by non- vanishing angular momentum 
J, which is obtained, for fixed /3, , by modifying the direction 
of the particles velocity vectors. The system total angular 
momentum is quantified by the spin parameter 

3/2||T|| 



A 



\Eu 



We consider A in the range 1 1.7 X 10"^ - 1.3 X 10"^ (see 
Table 1 for details). 



where Kt and 



^ We recall that typical values of A inferred from observations are 
~ 10"'^ for bright ellipticals up to ~ 10^^ for faint ellipticals (e.g. 
Lake 1983). 
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Table 1. Simulations parameters. 





R 


fit 






J V * 








A/10^^ 


c/a 


b/a 






TTlc 








pl04n0002 


0.002 


- 





- 


409600 


- 


P 


- 





0.47 


0.97 


4.3 ±0.06 


4.2 ±0.08 


3.6 ± 0.1 


0.52 


0.54 


0.05 


pl04n0005 


0.005 


- 





- 


409600 


- 


P 


- 





0.46 


0.93 


4.3 ±0.07 


4.2 ±0.1 


3.8 ± 0.1 


0.50 


0.53 


0.09 


plloOOl 


0.01 


- 





- 


102400 


- 


P 


- 


1.74 


0.44 


0.81 


4.6 ±0.08 


3.7±0.2 


3.6 ± 0.2 


0.43 


0.59 


0.26 


plSoOOl 


0.01 


- 





- 


102400 


- 


P 


- 


3.63 


0.55 


0.75 


5.2 ±0.2 


4.4 ±0.09 


3.7 ±0.2 


0.29 


0.49 


0.26 


pllOoOOl 


0.01 


- 





- 


102400 


- 


P 


- 


3.70 


0.56 


0.75 


5.0 ±0.2 


4.9 ±0.1 


3.9 ± 0.1 


0.27 


0.49 


0.28 


pl20o001 


0.01 


- 





- 


102400 


- 


P 


- 


3.75 


0.57 


0.80 


4.5 ±0.2 


4.4 ±0.1 


3.9 ± 0.1 


0.28 


0.47 


0.23 


pl20f001 


0.01 


- 





- 


102400 


- 


P 


- 


4.98 


0.24 


0.32 


5.1 ±0.2 


4.6 ±0.09 


3.9 ±0.09 


0.38 


0.53 


0.23 


pl04n001 


0.01 


- 





- 


409600 


- 


P 


- 





0.47 


0.97 


4.2 ±0.07 


4.1 ±0.09 


3.7 ±0.1 


0.50 


0.53 


0.05 


hq2r001 


0.01 


- 





- 


204800 


- 


1 


- 





0.47 


0.78 


5.7 ±0.3 


5.1 ±0.2 


4.7 ±0.2 


0.37 


0.51 


0.23 


pl20f002 


0.02 


- 





- 


102400 


- 


P 


- 


7.09 


0.49 


0.68 


4.6 ±0.1 


4.6 ±0.1 


4.0 ±0.06 


0.38 


0.53 


0.23 


pl2r002 


0.02 


- 





- 


204800 


- 


P 


- 





0.48 


0.87 


4.0 ±0.1 


3.8 ±0.1 


3.6 ± 0.1 


0.42 


0.51 


0.13 


hq2r002 


0.02 


- 





- 


204800 


- 


1 


- 





0.48 


0.88 


6.0 ±0.2 


5.7 ±0.3 


5.4 ±0.3 


0.41 


0.49 


0.12 


g02r002 


0.02 


- 





- 


204800 


- 





- 





0.49 


0.93 


6.2 ±0.2 


6.3 ±0.2 


5.3 ± 0.2 


0.49 


0.52 


0.10 


pl04n0025 


0.025 


- 





- 


409600 


- 


P 


- 





0.46 


0.63 


4.5 ±0.3 


3.8 ±0.1 


3.8 ± 0.2 


0.22 


0.50 


0.40 


pl20c005 


0.05 


- 





- 


102400 


- 


P 


- 


12.6 


0.49 


0.62 


4.7 ±0.1 


3.8 ±0.1 


4.0 ±0.05 


0.20 


0.53 


0.41 


pllnOOS 


0.05 


- 





- 


1024000 


- 


P 


- 





0.40 


0.66 


4.4 ±0.2 


4.0 ±0.1 


3.7 ±0.2 


0.35 


0.58 


0.37 


pl04n005 


0.05 


- 





- 


409600 


- 


P 


- 





0.41 


0.70 


4.3 ±0.2 


4.0 ±0.2 


3.6 ± 0.2 


0.38 


0.57 


0.33 


hq2r005 


0.05 


- 





- 


204800 


- 


1 


- 





0.41 


0.44 


8.1 ±0.5 


6.2 ±0.2 


6.2 ± 0.2 


0.11 


0.62 


0.58 


pl04n01 


0.1 


- 





- 


409600 


- 


P 


- 





0.96 


0.97 


4.0 ±0.3 


4.2 ±0.3 


3.8 ± 0.3 


0.01 


0.04 


0.03 


hq2r01 


0.1 


- 





- 


204800 


- 


1 


- 





0.39 


0.40 


7.4 ±0.6 


5.3 ±0.2 


5.2 ± 0.2 


0.02 


0.63 


0.64 


g02r01 


0.1 


- 





- 


204800 


- 





- 





0.42 


0.46 


7.5 ±0.5 


5.0 ±0.1 


5.2 ± 0.2 


0.08 


0.59 


0.56 


pl04n02 


0.2 


- 





- 


409600 


- 


p 


- 





0.96 


0.97 


3.6 ±0.3 


3.7 ±0.3 


4.0 ±0.4 


0.01 


0.01 


0.01 


hq2r02 


0.2 


- 





- 


204800 


- 


1 


- 





0.98 


0.99 


5.6 ±0.3 


5.6 ±0.3 


5.6 ± 0.3 


0.01 


0.02 


0.01 


hpl4r05m001 


0.01 


0.7 


0.5 


4 


163840 


81920 


p 


1 





0.28 


0.40 


4.3 ±0.2 


3.1±0.1 


3.1 ± 0.1 


0.35 


0.75 


0.63 


hpl4rlm001 


0.01 


0.8 


1 


4 


102400 


102400 


p 


1 





0.24 


0.34 


3.7±0.1 


2.8 ±0.07 


2.8 ±0.06 


0.31 


0.77 


0.68 


hpl20f4rlm001 


0.01 


0.8 


1 


4 


102400 


102400 


p 


1 


3.14 


0.24 


0.32 


3.8 ±0.1 


2.7 ±0.07 


2.7 ±0.06 


0.28 


0.79 


0.71 


g0pl04rlm001 


0.01 


0.3 


1 


0.4 


102400 


102400 


p 








0.37 


0.57 


3.3 ±0.05 


2.9 ±0.06 


2.9 ±0.05 


0.35 


0.64 


0.46 


gOplOSrlmOOl 


0.01 


0.3 


1 


0.8 


102400 


102400 


p 








0.43 


0.87 


4.3 ±0.09 


3.7 ±0.08 


3.6 ± 0.1 


0.48 


0.56 


0.13 


hpl025r2m002 


0.02 


0.4 


2 


0.25 


102400 


204800 


p 


1 





0.34 


0.35 


3.4 ±0.1 


3.1 ±0.03 


3.0 ±0.03 


0.0 


0.68 


0.68 


plpl05r2m002 


0.02 


0.1 


2 


0.5 


204800 


409600 


p 


p 





0.30 


0.34 


4.0 ±0.1 


3.4 ±0.04 


3.9 ± 0.1 


0.22 


0.71 


0.71 


plpl05r5m002 


0.02 


0.02 


5 


0.5 


81920 


409600 


p 


p 





0.47 


0.72 


4.8 ±0.1 


4.0 ±0.1 


4.8 ± 0.1 


0.32 


0.49 


0.29 


g0hq05r2m002 


0.02 


1 


2 


0.5 


204800 


204800 


1 








0.42 


0.85 


5.4 ±0.2 


5.6 ±0.3 


5.2 ± 0.2 


0.49 


0.57 


0.17 


hus20rlm002 


0.02 


0.9 


1 


20 


102400 


102400 


u 


1 





0.87 


0.90 


5.6 ±0.4 


6.5 ±0.6 


5.6 ± 0.4 


0.01 


0.11 


0.11 


hpl4rlm005 


0.05 


0.8 


1 


4 


102400 


102400 


p 


1 





0.98 


0.99 


2.5 ±0.03 


2.5 ±0.04 


2.5 ±0.04 


0.01 


0.01 


0.01 


hpl4r2m005 


0.05 


0.8 


2 


4 


102400 


204800 


p 


1 





0.99 


1.00 


2.2 ±0.03 


2.3 ±0.04 


2.2 ±0.03 


0.01 


0.00 


0.00 


hpl8r2m005 


0.05 


0.9 


2 


8 


102400 


204800 


p 


1 





0.99 


1.00 


2.0 ±0.03 


1.9 ±0.02 


1.9 ±0.02 


0.00 


0.01 


0.01 


hpl2r4m005 


0.05 


0.8 


4 


2 


81920 


327680 


p 


1 





0.99 


1.00 


2.1 ±0.03 


2.1 ±0.03 


2.1 ±0.03 


0.00 


0.01 


0.00 


npl4r2m005 


0.05 


0.9 


2 


4 


102400 


204800 


p 


N 





0.33 


0.33 


3.4 ± 0.2 


2.8 ± 0.09 


2.8 ± 0.08 


0.01 


0.01 


0.01 


npl4r4m005 


0.05 


1 


4 


4 


81920 


327680 


p 


N 





0.97 


0.98 


2.4 ±0.05 


2.3 ±0.04 


2.3 ±0.05 


0.01 


0.00 


0.00 


gl5pl4r2m005 


0.05 


1 


2 


4 


102400 


204800 


p 


1.5 





0.98 


0.99 


1.9 ±0.04 


2.0 ±0.05 


2.0 ±0.04 


0.01 


0.02 


0.00 


g0pl4rlm005 


0.05 


0.7 


1 


4 


102400 


102400 


p 








0.33 


0.33 


3.2 ±0.2 


2.8 ±0.06 


2.7 ±0.07 


0.01 


0.70 


0.70 


g0hq05rlm005 


0.05 


0.2 


1 


0.5 


102400 


102400 


1 








0.39 


0.39 


12.1 ± 1.2 


7.3 ±0.3 


7.8 ±0.4 


0.03 


0.64 


0.63 


hpl2r2m01 


0.1 


0.8 


2 


2 


102400 


204800 


p 


1 





0.98 


0.99 


2.6 ±0.1 


2.6 ±0.1 


2.6 ±0.1 


0.01 


0.03 


0.02 



First column; name of the simulation. /3* and /^i^: initial stellar and halo virial ratios, respectively. Mjj/M*: halo to stellar mass ratio. r«/rij; ratio of initial characteristic radii. Af» and A'^i^: 
number of stellar and halo particles, p* and pi^: initial stellar and halo distribution (P, N, U: Plummer, NFW, and uniform spheres distributions, respectively. Numbers identify 7 = 0, 1, 1.5 
models.). A: spin parameter, c/a and b/a: minor-to-major and intermediatc-to-major axis ratios, nia, rriij, rric and ea, ^c'- best-fit Sersic indices and ellipticities for projections along the 
principal axes. 
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2.1.2 Halo component 

In most of the simulations the dark matter halo component 
is represented by a 7-model distribution 

3-7 Mhrh 



Ph(r) = 



- 7 

47r r^{rh + r)* 



(8) 



where Mh is the total mass, rh is the characteristic radius, 
and we adopt 7 — 0, 1, or 1.5. In two simulations the 
halo component is represented by a Navarro, Frenk & White 
(1996, NFW) density profile with scale radius rh, truncation 
radius lOrh, and total mass Mh. Note that for the 7 = 1 
model (Hernquist 1990) ph{r) ~ for r 0, as in the 
NFW density profile. The velocity distribution of the halo 
particles is isotropic and, with few exceptions, such that the 
halo would be in equilibrium in the absence of the stellar com- 
ponent. In practice, in most cases the halo is nearly-virialized 
(0.7 < /3h ^ 1), while in a subset of simulations the halo 
is strongly out of equilibrium (/3h <^ 1; see Table 1 and Sec- 
tion 3.4). Length-scale ratios r,/ri^, and mass ratios Mi^/M, 
characterizing the initial conditions are listed in Table 1. 



2.2 Numerical methods 

Numerically, the initial distribution of particles in phase- 
space is realized as described in Nipoti, Londrillo & Ciotti 
(2002). For the stellar distribution we use typically ~ 
1 — 4 X 10"" particles, and in the two-component cases the 
number of particles A'h of the dark matter distribution is 
such that halo and stellar particles have the same mass (with 
few exceptions in which halo particles are twice as massive 
as stellar particles; see Table 1). 

For the simulations we used the parallel N-body code 
FVFPS (Fortran Version of a Fast Poisson Solver; Lon- 
drillo, Nipoti & Ciotti 2003; NLC03), based on Dehnen's 
(2000, 2002) scheme. We adopt the following values for the 
code parameters: minimum value of the opening parameter 
Omin = 0.5, and softening parameter e — 0.005 — 0.025 in 
units of r, (depending on the number of particles) . The time- 
step At, which is the same for all particles, is allowed to vary 
adaptively in time as a function of the maximum particle den- 
sity pmax. In particular, we adopted At = 0.3/ y' (47rGpmax), 
the classical time-step criterion for the stability of the leap- 
frog time integration. This gives time-steps as small as 
6 X 10~^tdyn during the first phases of the coldest collapses, 
assuring energy conservation to within 10~^ — lO"**. In order 
to estimate discreteness effects, we ran several test simula- 
tions varying the softening length e and the number of par- 
ticles. We found that the main properties of the end-product 
of a simulation with A'^, ~ 10® are not significantly differ- 
ent from those of lower-resolution simulations with the same 
initial conditions. In particular, the tests indicate that, for a 
one-component simulation with ~ 4x 10^ and £ = O.Qlr*, 
the final density distribution can be trusted down to radii 
~ 2e (~ 0.01 tm, where tm is the half- mass radius), in ac- 
cordance with numerical convergence studies (e.g. Power et 
al. 2003). 

The intrinsic and projected properties of the end- 
products are determined following Nipoti et al. (2002), while 
the position of the centre of the system is determined using 
the iterative technique described by Power et al. (2003). In 
particular, we measure the axis ratios c/a and b/a of the iner- 
tia ellipsoid (where a, b and c are the major, intermediate and 
minor axis) of the stellar component, their angle-averaged 
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Figure 1. Angle-averaged spatial density profile of the end- 
products of one-component collapses (panel a; solid lines) and of 
the stellar component of two-component collapses (panel b; solid 
lines). The dashed line is the 7 = 1.5 profile. Panels c and d: 
projected density profiles (along the principal axes) of the stellar 
component of the end-products of one- and two-component col- 
lapses. The dashed line represents the i?'^/* law. Different colours 
indicate different initial stellar density profiles: Plummer (blue), 
7 = 1 model (red), 7 = model (green), hus20rlm002 (magenta, 
see Table 1). Blue, red, green and magenta curves are vertically 
offset by —0.3, 0.3, 0.6 and 0.6 dex (top panels), and by 1, — 1, 
—2, and —2 magnitudes (bottom panels), respectively. 



density distribution and half-mass radius vm. For each end- 
product, in order to estimate the importance of projection 
effects, we consider three orthogonal projections along the 
principal axes of the inertia tensor, measuring the ellipticity 
e = 1 — 6e/oe, the circularized projected density profile and 
the circularized effective radius Rc = Vocbc (where Oc and be 
axe the major and minor semi-axis of the efi'ective isodensity 
ellipse). We fit (over the radial range 0.1 < R/Rc ^ 10) 
the circularized projected stellar density profile of the end- 
products with the ii^/™ Sersic (1968) law: 
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Figure 2. Circularized projected stellar density profiles of the 
end-products of representative one-component (left) and two- 
component (right) collapses. The dotted lines are the best-fitting 
Sersic models, l-cr error bars are also plotted. 



I{R) — Ic exp < —b(m 



(i?o) 



l/n 



- 1 



(9) 



where le = I{Rc) and b(m) ~ 2m- l/3-|-4/(405m) (Ciotti & 
Bertin 1999). In the fitting procedure m is the only free pa- 
rameter, because Tie and 7e are determined by their measured 
values obtained by particle count. 



3 PROPERTIES OF THE END-PRODUCTS 
3.1 Density distribution 

3.1.1 Spatial density profile 

In Fig. ^ we plot the angle-averaged density profiles of the 
end-products of one-component collapses (solid lines). The 
relatively small spread of the curves in the diagram indicates 
that the profiles are similar over most of the radial range 
(0.1 < r/ru 10), independently of the details of the ini- 
tial density distribution and of the specific angular momen- 
tum. In particular the profiles are reasonably well described 
by the 7 = 1.5 model (dashed line). 

In Fig. ^ we plot the angle-averaged density profiles 
of the stellar component of two-component simulations. The 
final stellar density profiles of two-component collapses are 
more varied than the results of one-component collapses, and 
they are typically flatter ai r < ru and steeper at r > tm 
than the 7 — 1.5 model. So the overall (and expected) effect 
of the presence of the halo in the collapse is to introduce 
significant structural homology. 

The initial and final halo density distributions are very 
similar when the halo is dominant {Mh/Mt > 2). When 
Mh/M, < 1 the halo profile in the central regions is signifi- 
cantly modified by the collapse. In particular, the final inner 
halo profile is shallower than the initial one if r,/rh > 1 
at t — 0, and steeper if Vt/rh ^ 1- Thus, the collapsing 
stellar component may be able to modify the density distri- 
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Figure 3. Left: best-fit Sersic parameter m as function of the 
initial dark-to-luminous mass ratio (''M ) /^* (''M ) within the 
half-mass radius of the stellar distribution. Right: best-fit Sersic 
index as a function of the ellipticity. For each model we plot three 
points corresponding to the three principal projections. Colours 
have the same meaning as in Fig.lTI 



bution of a pre-existing halo^. The final total (stellar plus 
halo) density profiles are characterized by significant struc- 
tural non-homology, and are typically well represented by 
7-models with < 7 < 2. 

Finally, we note that in both one- and two-component 
simulations the final density distribution tends to be steeper 
when the initial stellar distribution is a 7 = (green curves 
in Fig. or a 7 = 1 (red curves) model than for Plummer 
initial conditions (blue curves; see also Section 3.4). 



3.1.2 Projected density profile: Sersic fit 

The projected density profiles of the stellar end-products 
are analysed as described in Section 12.21 The fitting ra- 
dial range 0.1 < R/Rc 10 is comparable with or 
larger than the typical ranges spanned by observations (e.g. 
Bertin et al. 2002). The best-fit Sersic indices rUa, mt and 
rric (for projections along the axes a, b, and c, respectively) 
are reported in Table 1, together with the la uncertainties 
corresponding to Ax^ ~ 1. We note that the relative uncer- 
tainties on the best-fit Sersic indices are in all cases smaller 
than 10 per cent. 

Figure^ shows the circularized projected density pro- 
files of the end-products of one-component simulations (three 
projections for each system, each normalized to the corre- 
sponding value of 7e), while the de Vaucouleurs profile is 
the dashed straight line. Apparently all the end-products of 
one-component collapses do not deviate strongly from the 
R^^"^ law over most of the radial range, in agreement with 
previous studies (vA82, Londrillo et al. 1991, Trenti et al. 
2004). In particular, the Sersic index is found in the range 



^ Analytical restrictions on the relative density radial trend of 
stars and dark matter have been derived in Ciotti & Pelle- 
grini (1992), and Ciotti (1996; 1999). 
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Figure 4. Top: angle-averaged density profiles of the end-products 
of ono-component simulations with /3* = 0.002, 0.005, 0.01, 0.025, 
0.05, 0.1 and 0.2 (solid lines). The arrow indicates the typical 
simulation softening-length s, while dotted curves indicate l-cr 
uncertainties due to statistical noise. Bottom: circularized pro- 
jected density profiles of two representative one-component col- 
lapses. The bars are l-cr uncertainties. The dotted lines are the 
best-fitting Sersic models over the radial range 0.1 < R/Rc < 10. 



3.6 rn <j 8. The quality of the fits is apparent in Fig. |5| 
(left), which plots the surface brightness profile of a projec- 
tion of one of these end-products, together with the best-fit 
(m — 4.2 ± 0.07) Sersic law, and the corresponding residu- 
als. The average residuals between the data and the fits are 
typically 0.04 < (A/i) < 0.2, where ^l = -2.5 log /(i?)//o. 

In Fig. ^ we finally plot the circularized projected stel- 
lar density profiles of the end-products of two-component 
simulations. These systems deviate systematically from the 
7?^/* (straight line), and in most cases the profile remain be- 
low it at small and large radii. For these systems we found 
1.9 < m < 12, with average residuals in the same range as 
those of one-component collapses. We note that values of m 
significantly larger than 4 are found only for collapses start- 
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Figure 5. Break radius normalized to the effective radius as a 
function of the initial virial ratio for the simulations in Fig. |^ (top) 
with 0t < 0.1. The symbols refer to the average among the three 
considered projections values, which span the range represented 
by the vertical bars. 
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Figure 6. Final axis ratios c/a as a function of b/a (left), and of 
the initial stellar virial ratio /3« (right). The meaning of symbols 
and colours is the same as in Fig. |21 



ing from sufficiently concentrated stellar distributions (see 
Section 3.4). Figure |5| (right) plots the projected profile of a 
representative two-component simulation with Mh/M, — 2 
together with the best-fit (m — 2.2 ± 0.03) model and the 
residuals. 

The situation is summarized in the left panel of Fig. El 
where it is apparent how the initial amount and concentra- 
tion of the dark matter is crucial in determining the final 
shape of the profile. In particular, the trend is that m is 
anti-correlated with the initial ratio Mh{rM)/Mt{rM) within 
r-M- However, it is also apparent from Fig.|3]that there is sig- 
nificant spread in m for similar values of Mh(rM)/Af»(rM), 
which reflects differences in the initial stellar density profiles 
and halo virial ratios (see Section 3.4). We also note that for 
some strongly aspherical end-products rua, mt, and rric are 
significantly different (see Table 1). 
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Figure 7. Density, circular velocity, radial velocity dispersion, and anisotropy parameter for the end-products of a one-component sim- 
ulation (left panels) and a two-component simulation (right panels). Solid lines refer to stars; dashed lines to dark matter. Velocities are 
normalized to W GMt /tm- 



3.1.3 Behaviour at small radii 

So far we have considered the properties of intrinsic and pro- 
jected density profiles at radii > O.lrM (and O.liic)- We 
now focus on the behaviour of the profiles at smaller radii, 
yet confining our discussion to r > 0.01 tm where artificial 
numerical effects do not affect the profile (see Section f2. 211 . 

As apparent from Fig. |1| (top), the end-product density 
profiles of one-component collapses have fiat cores at r < 
O.lrM, in agreement with previous studies (vA82; May & 
van Albada 1984). Correspondingly, the projected density 
profiles are characterized by a break radius _Rb, in the sense 
that for R < Rh they stay below the best-fit Sersic profiles 
over 0.1 < R/Rc < 10, as shown in Fig. 2] (bottom panels), 
for a case with /?* = 0.002 (left) and a case with /3« = 0.02 
(right). Thus, /?* determines the radial range over which the 
final surface density profile is well fitted by the Sersic law, 
with colder initial conditions producing smaller cores (see 
Fig. ISJ. Note that we do not represent Rh/Rc for /9, = 0.2, 
because in this case the core extends to _R > 0.1-Rc, affecting 
the Sersic fit in the considered radial range 0.1 < R/Rc < 10, 
so the break radius is not well defined. 

Also the projected density profiles of two-component 
end-products fiatten at small radii and deviate from an in- 
wards extrapolation of the best-fitting Sersic law. However, 
the fiattening is typically more gradual and the break is not 



as apparent as in the one-component cases. This is in part due 
to the fact that the stellar end-products of two-component 
simulations are characterized in general by smaller m, corre- 
sponding to quite shallow profiles. 



3.1.4 Intrinsic shape and projected ellipticity 

Figure |S] shows the location of the end-products in the 
(6/a, c/a) and (/3t, c/a) planes. The final states of one- 
component collapses are roughly spherical for /3* > 0.1, pro- 
late (c/a ~ b/a ~ 0.3) for 0.05 < /?, < 0.1 and oblate 
{c/a ~ c/b ~ 0.5) for /3« < 0.05, in agreement with the find- 
ings of previous works (Aguilar & Merritt 1990; Londrillo 
et al. 1991). A similar behaviour is found in two-component 
simulations, though in this case the halo mass and concentra- 
tion, beside (3,, play a role in determining the final shape of 
the stellar component. In particular, in case of dominant halo 
the remnant is roughly spherical even for /3« = 0.02 — 0.05. 
When the stellar component develops a bar, the dark matter 
halo becomes mildly flattened (c/a ~ b/a ~ 0.7 — 0.8). 

The ellipticity e (measured for each of the principal pro- 
jections of each end-product) is found in the range — 0.8, 
consistent with those observed in real ellipticals, and no sig- 
nificant correlation with the Sersic index m is found (Fig. |3 
right panel). In particular, m > 4 systems span a wide range 
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Figure 8. Line-of-sight velocity (dotted line) and velocity dis- 
persion (solid line), along the principal axes, normalized to 
\J GMt I Rc , for the stellar component of four end-products with 
different values of A/h and A. 



in ellipticity, while lower-m systems are found both round 
and very flattened (though m as low as 2 is found only for 
e~0). 



3.2 Kinematics 

The internal dynamics of the end-products is quantified by 
a 'bona fide' circular velocity = G'Mtot(r)/r, the (angle- 
averaged) radial component a,, of the velocity dispersion ten- 
sor, and the anisotropy parameter a{r) = l — 0.5a^/a^, where 
at is the tangential component of the velocity-dispersion 
tensor. These quantities, together with the corresponding 
angle-averaged density distributions, are shown in Fig. Q 
for two representative simulations. The radial dependence 
of Vc , Or and a is qualitatively similar in the one- and two- 
component cases. The circular velocity is roughly constant 
for 0.2 < r/rM ^ 1 and decreases inwards and outwards; Ur 
has a plateau at small radii and decreases at larger radii. The 
systems are approximately isotropic in the centre [ct ~ 0) and 
strongly radially anisotropic for r > tm, in agreement with 
previous results (e.g. vA82; Trenti et al. 2004). 

For each model projection we compute the line-of-sight 
streaming velocity vios and the line-of-sight velocity disper- 
sion (Jias, considering particles in a slit of width -Ro/4 around 
the semi-major axis of the isophotal ellipse. Figure|H|plots vxos 
(solid line) and aios (dashed line) of the stellar component for 
four systems with different combinations of dark matter mass 
Mh and spin parameter A. The projected velocity dispersion 
profiles do not present the central depressions characterizing 
the isotropic Sersic profiles or the Hernquist model. In the 
presence of significant angular momentum (bottom panels) 
uios becomes comparable with ctios at _R ~ few Rc when the 
line-of-sight direction and the angular momentum vector are 
orthogonal. 



Figure 9. Left panels: initial (top) and final (bottom) differential 
energy distribution for a representative one-component simulation. 
Right panels: same as left panels for a two-component simulation. 
Solid, dashed and dotted curves refer to the stellar, halo, and total 
distributions, respectively. The energy per unit mass E is normal- 
ized to -Etot/Mtot. 

3.3 Differential energy distribution 

In Fig. El (left panels) we plot the initial (top) and final (bot- 
tom) differential energy distributions N{E) (i.e. the number 
of particles with energy per unit mass between E and E+dE) 
for a representative one-component collapse simulation. The 
final differential energy distribution is well represented by an 
exponential N{E) oc exp (tjE), over most of the energy range, 
with r; ~ 2.1 (when the energy per unit mass is normalized 
to GM,/Rc, see Binney 1982; van Albada 1982; Ciotti 1991). 
We find that the colder the initial conditions, the larger the 
range in E in which N{E) is exponential, because only very 
cold initial conditions produce strongly bound particles (and, 
correspondingly, populate the very central regions; see Sec- 
tion 13. 1.311 . 

In the right panels we plot, for a representative two- 
component collapse, the initial and final N{E) of the stel- 
lar (solid lines) and halo (dashed lines) component, together 
with the total N{E) (dotted lines). Note how the energy dis- 
tribution of the final stellar component is steeper than the 
halo component, i.e., how the very-bound end of the energy 
distribution (typically corresponding to particles in the cen- 
tral region of the density distribution) is populated by halo 
particles (dashed line), while most of the stellar particles are 
close to i? ~ 0. Interestingly, the slope of the total differen- 
tial energy distribution (dotted line) is very similar to that 
found in the one-component case when E is normalized to 
-Etot/Mtot. 

3.4 Exploring different initial conditions 

In general, the large-scale properties of the end-products 
of one-component simulations do not depend strongly on 
the initial distribution, provided that it is clumpy and cold 
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enough. However, some information on the initial distribu- 
tion is not completely erased during the collapse. As appar- 
ent from Fig. |3 while Plummer (blue symbols) initial con- 
ditions lead to systems very close to the de Vaucouleurs law 
(3.6 ^ m < 5.2), the end-products of 7 = (green symbols) 
and 7 = 1 (red symbols) initial distributions are character- 
ized by higher Sersic indices (5 < m < 8). 

In two-component simulations the initial shape of the 
dark matter halo does not have strong influence on the prop- 
erties of the stellar end-product. In particular, for a wide 
range of halo density profiles (7 = 0, 7 = 1,7= 1.5, and 
NFW distributions), the final stellar density profiles have 
m < 4, down to 2, if the halo is sufficiently concentrated. If 
the halo is light and/or extended, independently of its shape, 
the final stellar components have m > 4. We also consid- 
ered cases in which the halo, as well as the stellar compo- 
nent, is strongly out of equilibrium (/3h ^ 1; see Table 1). 
Independently of the details of their initial halo density dis- 
tributions, the stellar end-products of these simulations have 
global properties more similar to one-component collapses 
than to collapses in nearly virialized haloes. As a conse- 
quence, these end-products (represented as stars in Fig. |3J, 
contribute significantly to the scatter in the left diagram of 
Fig. El because they typically have higher values of m than 
the collapses in nearly virialized haloes with similar dark-to- 
luminous mass ratio. In addition, in analogy with the one- 
component case, two-component collapses with 7 = 1 initial 
stellar distribution (such as g0hq05r2m002, red filled circles 
in Fig. I^J produce systems with higher Sersic index than 
those with Plummer initial distribution (blue symbols) with 
similar dark-to- luminous mass ratio (see Fig. |3 left panel). 
Interestingly, the stellar system with the highest Sersic in- 
dex [m ~ 12.1) was obtained starting with a 7 = 1 stellar 
distribution embedded in a collapsing dark matter halo (sim- 
ulation g0hq05rlm005, red stars in Fig.l^J. 

We also considered a rather extreme case of two- 
component initial conditions (model hus20rlm002 in Table 1, 
magenta curves and symbols in Fig. 13| and [SJ that could 
be also interpreted as a multiple merging experiment. The 
halo is a Hernquist model with core radius rh, and the stars 
are distributed in 20 uniform spheres of radius 2rh, whose 
centres of mass are uniformly distributed in a sphere of ra- 
dius 20rh. Each small sphere is cold, having all particles at 
rest with respect to its centre of mass. So the small spheres 
collapse first because of self-gravity and then merge because 
of dynamical friction, finally forming a stellar system with 
best-fit Sersic index^ m ~ 6.5. We note that in this case the 
final central halo density distribution is significantly fiatter 
than the the initial one, as expected as a consequence of dy- 
namical friction heating (see also Nipoti et al. 2004). 

According to these results, it appears that the only way 
to produce m > 5 systems with dissipationless processes is 
either merging or collapse starting from concentrated and 
cold initial distributions, while dissipationless collapses of 
more diffuse, cold distributions lead to m < 5 systems. We 
note that lower values of m were also found as a result of 
quite different dissipationless processes: Merritt et al. (2005) 
found that the projected density of the dark matter haloes of 
cosmological N-body simulations are well fitted by the Sersic 



^ Similar values of m were also found by Nipoti et al. (2003b) for 
the end-products of galactic cannibalism simulations. 




Figure 10. Dark-to-total mass inside r for the end-products of 
the two-component simulations. Dotted, solid, and dashed curves 
refer to systems with average best-fitting Sersic index (m) < 3, 
3 < (m) < 5, and (m) > 5, respectively [(m) = (ma-|-m(,+mc)/3]. 



law with 2 < m < 4, while Nipoti et al. (2002) obtained 
1 5 for the end-products of radial orbit instability. 



4 DISCUSSION AND CONCLUSIONS 

Before summarizing the main results of the numerical simu- 
lations, it is useful to list three (of the several) observational 
properties of ellipticals that should be explained by any plau- 
sible formation (or evolutionary) scenario: 

(i) The best-fit Sersic index increases from m ~ 2 to m ~ 
10 from faint to luminous galaxies (Caon et al. 1993; Berlin 
et al. 2002). Thus, the light profiles of luminous ellipticals as 
a function of R/Re are steeper for R < Rc and flatter for 
R > Re than the profiles of low-luminosity ellipticals (see, 
e.g., Ciotti 1991). 

(ii) However, the light profiles of luminous ellipticals 
present a flat core in their very central regions, deviating from 
the global best-fit Sersic law. The profiles of low-luminosity 
ellipticals are instead well fitted by the Sersic law down to 
the centre (e.g. Lauer et al. 1995; Trujillo et al. 2004), so 
in the very central regions luminous ellipticals have fiatter 
density distributions than low-luminosity ellipticals . 

(iii) Finally, luminous ellipticals have lower mean density 
within Rc than low- luminosity ellipticals, a well known con- 
sequence of the Kormendy (1977) relation. 

In this paper we presented the results of numerical N- 
body simulations of one- and two-component dissipationless 
collapses, focusing on their possible relevance to the origin 
of the weak homology and of the core/power-law dichotomy 
of ellipticals, i.e., we addressed points (i) and (ii) above, 
while it is quite obvious that point (iii) cannot be studied 
with dissipationless collapse simulations. Overall, our results 
suggest that dissipationless collapse is effective in producing 
weak structural homology, kinematic properties consistent 
with observations, and cores on scales remarkably similar to 
those observed. In particular: 
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• The end-products of one-component dissipationless col- 
lapses typically have projected surface brightness profile close 
to the de Vaucouleurs model (see also vA82; Londrillo et 
al. 1991; Dantas et al. 2002; Trenti et al. 2004). When fitted 
with the Sersic law over the radial range 0.1 < R/Rc < 10, 
the resulting profiles are characterized by index 3.6 < m < 
8; m > 5 final states are obtained only for rather concen- 
trated initial conditions. 

• The end-products of collapses inside a dark matter halo 
present significant structural non-homology. The best-fit Ser- 
sic indices of the stellar projected surface density profile span 
the range 1.9 < m < 12. Remarkably, the parameter m cor- 
relates with the amount of dark matter present within Re, 
being smaller for larger dark-to- visible mass ratios. 

• The projected stellar density profiles are characterized 
by a break radius 0.01 < Rh/Re ^ 0.1 within which the pro- 
file is fiatter than the inner extrapolation of the global best- 
fit Sersic law. Colder initial conditions lead to end-products 
with smaller 7?b/^?e; in general, the resulting 'cores' are bet- 
ter detectable in high-m systems. 

The results above indicate that dissipationless collapse 
is able to produce systems with projected density profiles 
remarkably similar to the observed surface brightness pro- 
files, with high quality one-parameter Sersic fits even for 
low m values when non-negligible amounts of dark mat- 
ter are present. In particular, in Fig. IIUI we plot the (inte- 
grated) dark-to-total mass ratio Mh{r) /Mtot{r) for the two- 
component simulations end-products. A clear dichotomy is 
apparent: while in low-m systems (dotted lines) the mass- 
to-light ratio increases significantly inwards, the opposite is 
true for systems with m > 3 (with few exceptions). We also 
note that, while for m > 3 systems the dark matter mass 
inside the half mass radius is of the same order as the vis- 
ible mass, in accordance with observations (e.g. Berlin et 
al. 1994; Gerhard et al. 2001; Magorrian & Ballantyne 2001; 
Romanowsky et al. 2003; Treu & Koopmans 2004; Cappellari 
et al. 2006; Dekel et al. 2005; Samurovic & Danziger 2005), 
for very low-m systems the dark matter is dominant. 

According to the standard interpretation, the central 
cores observed in several bright ellipticals are a consequence 
of formation through merging, being produced by the inter- 
action of binary supermassive black holes with a stellar cusp 
(Makino & Ebisuzaki 1996; Faber et al. 1997; Milosavljevic 
& Merritt 2001). Our results, in agreement with previous 
numerical explorations (vA82; May & van Albada 1984), in- 
dicate that a break in the profile at small radii and a flat 
central core are features produced naturally by dissipation- 
less collapse. In particular, the fact that the size of the core is 
correlated (and the maximum central density anti-correlated) 
with the initial virial ratio is a direct consequence of Liouville 
Theorem* (e.g.. May & van Albada 1984). The higher resolu- 
tion of our simulations allowed us to investigate this effect in 
regions (down to i? ~ 0.01-Rc) comparable to those explored 
by high-resolution photometry of real ellipticals. In particu- 
lar, the profiles of the end-products of our single-component 
simulations are remarkably similar to those of observed 'core' 

* Hozumi, Burkert & Fujiwara (2000) pointed out that another 
consequence of the conservation of the phase-space density is that 
the initial anisotropy in the velocity distribution affects the fi- 
nal density profile at small radii. As shown by their spherically- 
symmetric dissipationless collapse simulations, smaller central 
cores are produced by more radially anisotropic collapses. 



ellipticals. For instance, it is suggestive to compare the pro- 
file plotted in Fig. |1] (bottom, left panel) with that of the 
core-elliptical NGC3348 (see fig. 10 of Graham et al. 2003), 
for which Trujillo et al. (2004) report best-fit m ~ 3.8 and 
Rh/Rc — 0.016. The fact that the end-products of our simu- 
lations reproduce very nicely the observed cores makes dissi- 
pationless collapse a plausible alternative to the binary black 
hole scenario for the origin of the cores. The next step of the 
present study would be the exploration with a hydro N-body 
code of the first stages of galaxy formation dominated by 
gaseous dissipation and the following feedback effects due to 
star formation. 
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